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Abstract 



4 supersymmetric theories. We present some = 1,2 Yang-Mills and Wess-Zumino models, 

m : 

■ in which the planar one-loop anomalous dimension matrix in the sector of holomorphic scalars 

' is identified with the Hamiltonian of an integrable quantum spin chain with SU (3) or SU{2) 

symmetry, even if the theory is away from the conformal points. This points to a more universal 
origin of the integrable structure beyond superconformal symmetry. We also emphasize the role 
^ , of the superpotential in the appearance of the integrable structure. The computations of operator 

I mixing in our examples by solving Bethe Ansatz equations show some new features absent in = 4 

SYM. 
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I. INTRODUCTION 



Whether the integrable structure that is abundant in two dimensional field theories or 
statistical models could emerge in four dimensions has been a fascinating topic on the re- 
search frontier. Indeed, a (non-topological) Yang-Mills theory in four dimensions should 
have very rich physics, so that it is hard to believe it could be exactly solved as a whole 
the same way as in two dimensions. On the other hand, one cannot either rule out the 
possibility that an integrable structure appears in a subsector or in a special limit. Actually, 
there have been quite a few evidences for an integrable structure in the self-dual sector of a 
Yang-Mills theory Q, fl J fl fl 3 E | or in the high energy limit of QCD 3, 0, Q- 

Recently there have been revived interests in the integrable structure in A/" = 4 super- 
conformal Yang- Mills (SCYM) theory in four dimensions, at least in the planar limit. One 



development 



ll\ was to apply the integrable structure, originally found in the QCD context 



12j |. to the computation of anomalous dimensions of twist-two operators in A/" = 4 SCYM. 
(This result was also obtained by means of other methods Q].) Another strong evidence 
came from a recent one-loop calculation of the anomalous dimensions of composite oper- 
ators 3|- (The interest in composite operators in the gauge theory was mainly inspired 
by the proposal Q] that at least some of them could be viewed as the holographic dual 
of string states in a special limit of the curved background AdS^ x S^.) It turned out 
that the planar one-loop mixing matrix for the composite operators, that consist of a string 
of scalar fields in the theory, is equivalent to the Hamiltonian of an integrable spin chain 
with SO {6) symmetry! One may naturally ask whether a similar integrable structure could 
appear in other Yang-Mills theories, or Wess-Zumino models, with less supersymmetries or 
even with conformal symmetry broken. This is an interesting issue, because one wants to 
know whether the appearance of the integrable structure in the A/" = 4 SCYM is related to 
the maximal superconformal symmetry that the theory has. 

In this paper, we present examples in which a similar integrable structure survives de- 
formation of the theory with less unbroken supersymmetries. More concretely, we show 
that in some A/" = 1 (and Af = 2) supersymmetric Yang-Mills theories in four dimensions, 
the planar one-loop anomalous dimension matrix for composite operators of holomorphic 
scalars is equivalent to the Hamiltonian of an integrable quantum spin chain, even if the 
theory is away from the conformal points. If we take a limit in which the 't Hooft gauge 
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couplings vanish, these models reduce to Wess-Zumino models, and the integrable structure 
still survives this limit. Our results indicate that the appearance of the integrable structure 
in Yang-Mills theory in four dimensions should have a more universal and profound origin, 
not restricted only to theories with (maximal) superconformal symmetry. Moreover, in our 
examples, the superpotential term is seen to play a crucial role in determining the integrable 
structure. 

Some of our models are deformation of orbifolding daughters of the A/" = 4 SCYM with 
the superpotential strength changed. In their conformal phase, they are reduced to a quiver 
theory with a product gauge group obtained by orbifolding the A/" = 4 SCYM with a 
discrete subgroup of SO (6) i?-symmetry. We therefore expect that there is a close relation- 
ship between integrable spin chains in these A/" = 1, 2 SYM and the 5*0(6) chain appearing 
in the Af = 4 SCYM. Indeed it turns out that all the integrable spin chains revealed by us in 
these A/" = 1, 2 SYM is closely related to a closed subsector of an SO (6) chain; at conformal 
points the parameters of the latter reduce to those of the 5*0(6) chain in A/" = 4 SCYM. 
However, the operator mixing that results from diagonalizing the spin Halmitonian is com- 
pletely different, because of an interesting interplay between the global symmetry index and 
the discrete index for gauge group factors, which is absent in the A/" = 4 SYM. 

The results for our A/" = 1, 2 models, combined together, provide a description of cascade 
breaking of the symmetry of the integrable spin chains starting from the A/" = 4 SYM. This 
motivates to make the conjecture that in all A/" = 1, 2 orbifolded daughters [3, ll^ of the 
A/" = 4 SYM (with a non-abelian global symmetry), with their 't Hooft and superpotential 
couplings deformed away from the conformal points while keeping the global symmetry 
of the superpotential, there is always a non-trivial integrable structure in operator mixing 
at planar one-loop level for (anti-) holomorphic composite operators, at least for those 
consisting of purely scalars without derivatives. Moreover, this integrable structure survives 
in the resulting Wess-Zumino models when all 't Hooft couplings are sent to zero. 

The contents of this paper are as follows. In section 2 we construct an A/" = 1 SYM model 
whose conformal phase is the orbifolding limit of A/" = 4 SCYM. In section 3 we compute the 
planar one-loop corrections to composite operators consisting of scalar fields of the model. 
In section 4 we obtain the mixing matrix for the renormalized composite operators, and 
show that the anomalous dimension matrix in the (anti-) holomorphic sector is equivalent 
to a Hamiltonian of an integrable SU{3) spin chain. We generalize our discussion to other 
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M = 1 (and M = 2) Yang-Mills and Wess-Zumino models in section 5, and devote section 6 
to a brief summary. In two appendixes we present examples for computing planar one-loop 
anomalous dimensions via solving the Bethe ansatz equations of the SU{2i) quantum spin 
chain in section 4 for our M = 1 model. This computation will explicitly demonstrate how 
and why the operator mixing differs from the Af = 4 SYM case, though the quantum spin 
chain can be viewed as a closed subsector of the latter. 



II. THE 7^ = 1 MODEL 



Various A/" = 1 superconformal gauge theories have been constructed 
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via consid- 



ering D3-branes at orbifold singularities of the form C^/F (with F a discrete subgroup of the 



S0{6) R-symmetry of A/" = 4 SYM) (l6|,|20|. In AdS/CFT correspondence j2l|, the F action 
is translated to an action AdS^ x (S^/T) with the AdS part unaffected. So the world- volume 
theory on D3-branes remains to be a conformal field theory. With appropriately chosen 
F, supersymmetry is broken down to A/" = 2, 1 or 0. 

The simplest example for A/" = 1 SCYM is the case with F = Z3 proposed in ref. jisl ]. 
The gauge group of the resulting theory has a product group U{N)^^^ x U{N)^'^^ x U{N)^^\ 
We will use a dsicrete index A = 1,2,3 to label these U{N) groups. The matter in the 
theory consists of the bi-fundamental chiral superfields 

3{{N, N, 1) © (iV, 1, N) © (1, N, N)}, (1) 

and their anti-chiral partners (Fig. 1). Inside each pair of circular brackets, we have the 
representations of the gauge groups, while the overall factor 3 in Eq. ([Q) reflects an SU{3) 
global symmetry inherited from the SO (6) symmetry of the Af = 4 SYM, that is broken 
down to SU{3) by orbifolding. 

Therefore, each of the U{N) groups is coupled to Nj = 3N fundamental chiral and 
Nf = 3N anti-fundamental anti-chiral matter. One may denote the bi-fundamental chiral 
matter fields in Eq. as ^"-'^^^^ Here the pair of indices, {A,B), labels which two of 
the gauge groups that the field is coupled to; and the index a = 1,2,3 is that for the 
representation 3 of SU{3). To simplify the notation, we will write $^ = <|)"'('^'^) with 
{A,B,C) being a cyclic permutation of (1,2,3); and denote their Hermitian conjugate as 
= (<l>^)"l' (Fig. 1), with the tilde symbol suppressed when this is no confusion. 
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FIG. 1: A quiver diagram for D3-branes on a C^/Zg orbifold. At the nodes 
we have vectormultiplets in the gauge group indicated, while the arrows 
connecting each pair of nodes correspond to the bi-fundamental fields. 

With the matter content (Q), one can construct an A/" = 1 gauge theory, which is not 
necessarily conformal invariant by allowing unequal gauge couplings and/or a more general 
superpotential. In the standard A/" = 1 superfield formalism, the Lagrangian of our model 
reads 



^ A 



1 

4 

i / £9ti{WAWA + h.c.) + ^ / d^e{W + h.c.) 



(2) 



In the first term, we assumed that {A,B,C) is a cyclic permutation of (1,2,3). Here to 
deform the theory obtained by orbifolding the A/" = 4 SYM, we allow the gauge couplings, 
qa with A = 1,2, 3, of the three U{N) groups to be different from each other, and take the 
superpotential W to be a generic SU{3) invariant one: 



h 



W = -e,fc,tr<l>?$^$^. 



(3) 



with h arbitrary. So the global symmetry remains to be SU (3) after the deformation. If the 
three gauge couplings are the same, our model also has a discrete Z3 symmetry acting on 
the index A. The one-loop beta functions in this model can be extracted from the NSVZ 

d\A 



I3h 



dlnfi 

dXh 
dlnfi 



^A j^]^ _ 



47r 



87r2 



7a, 



(4) 



where Xa = g\N {A = 1, 2, 3) are 't Hooft couplings for the three U{N) groups and 
Xh = \h\'^N; '-fA are anomalous dimensions of the three types of bi-fundamental scalar fields. 
In large N limit, the one-loop 7^ are 

^(E^^-2A.), (5) 



7a 



2 din a ' IGtt^ 



In this paper we will work in the region with weak 't Hooft couplings and small Xh- 
Substituting Eq. © into Eq. (jH) we have 

^ B 

''^-il^-^EA.-aM.otA'). (6) 

There is a line of fixed points, Ai = A2 = A3 = A/j, corresponding to the orbifolded SCYM. 
Though the above equations cannot be exactly solved we can consider the simplest case 
with Ai = A2 = A3, denoted as A, to demonstrate some features for the flow. First from 
the flow equations it is easy to see that if initially the three gauge couplings coincide, then 
they will keep to do so during the renormalization group flow. Secondly we can easily see 
that during the flow one always has Xh = cexp (Gvr^/A), with an arbitrary positive integral 
constant c. So for fixed value of c, both the 't Hooft and the superpotential couplings, A 
and Xh, respectively, possess a fixed point X = Xh = A*(c), whose value is determined by 
the root of the transcendental equation x = cexp (Gtt^/x). At the fixed point, the theory 
is an interacting conformal field theory at the quantum level. (The theory is no longer 
asymptotically free for finite A*(c)). When c varies one obtains a line of RG fixed points. 
When away from the conformal points, at infra-red the theory will be dominant by gauge 
interactions if A > A*, A/j < A* and by superpotential interactions if A < A*, Xh > X^. In 
the following we will calculate the operator mixing for small values of the couplings A^ and 
Xh. 



III. MATRIX OF RENORMALIZATION FACTORS 

We will study one-loop renormalization of composite operators, which are a product 
of bi-fundamental scalars 0^, the scalar component of the supermultiplets and their 
Hermitian conjugate 0^ = {(PaV, without derivatives: 

O^i-^i =tr (7) 
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where each index Ii {1 < I < L) stands for a pair of indices {ai,Ai) or ai,Ai. The index 
Ii for a fixed I can take 18 different values. The fields at the right-hand side are all taken 
to be at the same space-time point. These operators form (reducible) SU{3) tensors with 
L indices. In particular, this class of composite operators includes purely chiral (or holo- 
morphic) operators which are a product of 0^'s only (with no indices of type a). In physics 
one is restricted to composite operators that are gauge invariant after taking the trace over 
the gauge-group matrix indices. To compute operator mixing, it is better to work with 
the composite operators before taking the trace. We note that gauge invariance strongly 
constrain possible choice of the index sequence (/i,---,/^). (See below.) In general, the 
scalar operators ((Tj) mix under renormalization, and the renormalized operators with definite 
dimension are linear combinations of the bare operators. If we specify a particular operator 
basis, with / being a sequence /i,- ■■,//,, the matrix Zj of renormalization factors in 
this basis, 

OL = Z'jO', (/=!,. ..,18^), (8) 

is defined by the requirement of canceling UV divergences in the correlation functions 

{zI{'^Mx,)...zI/%,{xl) Ort-"^(x)), (9) 

These Zj depend on the UV cutoff A as well as various couplings of the theory (in the 
large- limit), and form a 18^ x 18^ matrix. In rest of this section, we will compute the 
6Zj = Z\ — 5 J due to planar one- loop diagrams, using the component field formalism as 
in ref. There are three types of planar one-loop diagrams, as shown in Fig. 2, that 

contribute to the correlations 0. We will choose the Fermi- Feynman gauge for gauge boson 
propagators, in which the anomalous dimension of a single scalar field is already given by 
Eq. (0). 

The one-loop self-energy diagram. Fig. 2c, leads to the wave function renormalization of 
scalar fields, which can be directly read off from Eq. (jSj). Half of the self-energy correction 
in correlation functions are cancelled by wave function renormalization of the external 
legs. The counterterms that cancel the remaining divergences are given by 

5Z^-:'l;^: = -1(7^, + 7a,,J in a S^:sf-l, (10) 

where each Ii or Ji stands for a pair of indices {ii,Ai) or {ji,Bi), respectively. Note that 
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a b c 

FIG. 2: One-loop planar diagrams on correlation functions (jH)). 

the contribution of this diagram to the Z-matrix is diagonal in the indices Ji and in 
h+i, Ji+i- 

As for the contributions from Fig. 2a and 2b, we have to distinguish between two different 
cases: i) two nearest-neighbor scalar fields both are holomorphic (or anti-holomorphic), 
namely (9 ~ ■ ■ ■ ' ' ' i'^^ ' ' ' ^'ai^'a^^i ' ' them is holomorphic and the 

other anti-holomorphic, i.e. (9 ~ ■ ■ ■ 'P'ai'P'aI'^ ' ' ' ' ' ' ' ' nearest-neighbor 

pairs in a composite operator O all belong to the above case i), we call it a holomorphic (or 
anti-holomorphic, respectively) operator. Otherwise, if the above case ii) happens to one 
nearest-neighbor pair, the composite operator is called non-holomorphic. 

The correction due to gauge boson exchange. Fig. 2a, contributes to the Z-matrix asso- 
ciated with holomorphic operators: 



5Z 



^ In A 5/' 5^+1 



-III, 



i-'i+i'" 



IGvr 



(11) 

where Af = (Ai ± 1) mod 3. For anti-holomorphic operators, one just replaces Aj^ A'^. 
Meanwhile, the same one-loop diagram. Fig. 2a, yields the following Z-matrix for non- 
holomorphic operators: 

A.- 



6Z 



{a.)-JiJi 



+1' 



-1,1, 



— ^lnA5/'5?+\ 



6Z 



{a)---JiJi^ 



A 



At 



h -'i+i 



-III, 



iJ-i+\--- 



(12) 



Note that here an overall factor of "^^,^,+1 is implied due to the requirement of 

gauge invariance. 

As for Fig. 2b, one has to be careful in extracting the SU{?)) structure of the resulting 
contributions to the Z-matrix, because the quartic scalar vertex in this diagram involves 
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both gauge and superpotential couplings. For holomorphic operators, we have 

^ _lnA I V _ ^ J A;,5/'+M/' }> . (13) 

For non-holomorphic operators the Z-matrix exhibits some comphcations: 

■■■/i/i+r-- g^2 I " h Ii+i ^ " 2 ^i'-'i+i 

In next section, we will consider the operator mixing resulting from above Z-factors. 



IV. OPERATOR MIXING AND SPIN CHAIN 

The anomalous dimension matrix (ADM) at one loop for operators (jZj) is determined by 
the standard arguments through 

r = (15) 

d\nA ^ ' 

Operator mixing arises when one diagonalizes this matrix to obtain operators with definite 
dimension. The ADM acts on a IS'^-dimensional vector space Vi ® ... ® Vl, where each Vi 
is a complex vector space spanned by 0^ and 0^, As proposed by Minahan and Zarembo 
jl^ for the A/" = 4 SYM, it is extremely useful to identify the ADM as the Hamiltonian of a 
lattice spin system, where the lattice sites are labeled by the subscript / = 1, 2, ■ ■ ■ , L. The 
main goal of our paper is to look for an integrable spin chain for the ADM, at least in a 
sector, in our model. 

However, there is a technical complication here. In the A/" = 4 case, the adjoint scalars 
are labeled only by one index i = 1, ■ ■ ■ , 6 for the 5*0(6) symmetry. But in our model, we 
need a pair of indices (a. A) (or a, A) to label the bi-fundamental scalars, where the extra 
index A is necessary to label the arrows in the quiver diagram, indicating which pair of 
U{N) gauge groups the scalar is coupled to. This is because our model is a deformation of 
an orbifolded model. We note, however, that in accordance with the quiver diagram. Fig. 



1, it is easy to see that gauge invariance of a composite operator (|7j) strongly constrains its 
index sequence Ji, ■ ■ ■ , 1^. For example, if J/ is known, then gauge invariance dictates the 
gauge index in J^+i, depending on whether its SU{3) index is of type a or type a. Thus if Ai 
is known, then gauge index sequence of a gauge invariant composite operator is completely 
determined by its SU{3) index sequence, corresponding to an SU{3) spin chain, with spin 
on each site belongs to either 3 or 3. The Hilbert space at each site is a 6-dimensional 
(reducible) representation space V of SU{3): = 3 © 3, and the dimensionality of the 
Hilbert space of the spin chain is thus reduced to 6^. 

In order to write ADM of operators (jZj) in compact form, and to facilitate later comparison 
with the Af = 4 mother SYM theory, we introduce the projection operators J^, which 
projects V to its invariant sub-spaces 3 and 3, respectively; in component form, they projects 
a vector f * = {v"', v"-) to its components v"' and to v"'. We also define the permutation 
operator P and the trace operator K, which act on the tensor product V ^V, respectively, 
as 

P{u ^ v) = V u, 

K{uiS)v) = {u ■ v)'^{ea<S) ea + ea<S) Co), (16) 

a 

where u, v are vectors in V and {u ■ v) = '^^i'^"''^'^ + u°'v°'); moreover, ia and 65 are vectors 
of an orthonormal basis in V. 

Then we write the ADM defined in Eq. ()15p as follows: 

L 

1=1 i,j=± 

n/+Mi) = ^lUAi) = + ^(1 - Pm+i), (17) 

r+;,(A) = + + Ki,i+i - K,,i-,iPi,i+i) - ^ - ^M+i)^M+i, 

1=1 1=1 

riUAi) = -jA, + + KiHi - Ki,i+iPi,i+i) - ^ - Ki,i+i)Pi,i+i. 

1=1 1=1 

We note that though the coefficients on the right side depend on the gauge indices A/, but 

the operators Pi^i+i and act only on the SU (3) indices. 

Since at planar one-loop level holomorphic (or anti-holomorphic) operators, consisting 

of 3's (or of 3's, respectively), mix only among themselves. So it makes sense to restrict 

10 



the ADM to the holomorphic (or anti-holomorphic) sector, and then to ask whether the 
so-restricted ADM can be identified with the Hamiltonian of an integrable spin chain or 
not. (From the quiver diagram, one can easily see that gauge invariance requires (anti- 
)holomorphic operators have length L = 3k with integer k.) 

From above equations we see that the ADM Fh, restricted in the holomorphic sector and 
associated with the projection operators J^" ■ ■ ■ J^, is given by 

L \ ^ 

1=1 1=1 

A same expression can be obtained for anti-holomorphic operators. 

It is worthy to note that in Eq. (fTHj) . the dependence on gauge indices Ai can actually 
be eliminated, since gauge invariance requires that L = 3k and the gauge index sequence 
goes around the triangular quiver diagram. Fig. 1, only in one direction. Therefore we can 
rewrite the ADM in the holomorphic sector as 



r(h) = To + Fi, 



L ' 



24vr2 



Xa — 3Xh 



A=l 

Ti = ^E(l-n..i)- (19) 
1=1 

It is easy to see that Fq is a constant, while Fi depends only on A^, the superpotential 
coupling. By introducing the spin operators 

S:^ = ^{St6'^'S]6l), (20) 

for each lattice site the ADM, Fi, can be rewritten in terms of manifest spin-spin interactions: 

L 

87r2" 



= ^ - iSfS^l,)'}. (21) 

Since only the permutation operator P appears in the ADM (|19p. it can be identified 
to be the Hamiltonian of an integrable spin chain with SU{3) symmetry, for arbitrary 't 
Hooft couplings A^. The detail on the integrability of SU{M) spin chain will be presented 
in appendix A. The Hilbert space of this spin chain is the tensor product Ui ® ■ ■ ■ ® Ul 
with Ui = spanned by covariant vectors of SU{3). Thus, though the full ADM (|T7j) in 
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our model does not in general correspond to an integrable spin chain, the ADM ()19p in the 
holomorphic sector does. 

It is interesting to note that in our model, an integrable SU (3) spin chain (fT^ appears 
even if the theory is away from the conformal fixed line, independent of the values of 
and A/j. The integrability of the system enable us to find the exact one- loop anomalous 
dimensions of (anti-) holomorphic operators via applying Bethe ansatz equations. Some 
examples and computational details are shown in appendix B. 

In Sec. II we have shown that at IR the gauge interaction and superpotential interaction 
are decoupled. In the limit A/i — 0, the ADM (fTn|) is proportional to identity operator so that 
(anti-)chiral operators ((Tj) do not mix with each other under planar one-loop renormalization. 
However, the case with A^ — > is more interesting. When A^ — > 0, the present M = 1 
gauge theory actually approaches to an A/" = 1 Wess-Zumino model. Then, according to 
the above Eq. ()19p. the planar one- loop ADM for composite operators consisting of purely 
(anti-) chiral scalars in this Wess-Zumino model also corresponds to an integrable SU{'i) 
spin chain. (Since the chiral supermultiplets are taken to be bi-fundamental representations 
of three global U (N) groups, here the planar limit makes sense as the limit in which N oo 
with hN'^ kept fixed.) We note that in this Wess-Zumino model there is no gauge invariance 
constraint, so the length L of the holomorphic composite operators for the integrable ADM 
(or quantum spin chain) do not need to be restricted to a multiple of three! 

To conclude, we make the following remark: Though in the above we have argued that 
the SU{3) spin chain arising from the ADM in the holomorphic sector is "blind" to the 
discrete index A, the operator mixing is not. Because each 0-field factor in the composite 
operators now carries the index A due to orbifolding, as we will see in appendix B, the 
operator mixing becomes completely different from that for the 5*0(6) spin chain in A/" = 4 
SYM. 



V. GENERALIZATION 
A. Lift to the 7^ = 4 SYM 

We have mentioned that on the conformal fixed line, our present A/" = 1 model is a Z3 
orbifolding of AT = 4 SYM Q Q. Then all correlation functions of the orbifolded theory 
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are known to coincide with those of TV = 4 SYM, modulo a rescahng of the gauge couphng 



constant. This has been shown either by using string theory 



24| . or by using Feynman 



diagrams in field theory 25|. Thus we expect that the spin Hamiltonian (|19j) is related to 



o y. 



the integrable 5*0(6) spin Hamiltonian obtained by Minahan and Zarembo 

To see this without explicit calculation, we may argue as follows. First notice that for 
fixed index A, the bi-fundamental 0" fields, together with their conjugate (pa, form a 6- 
dimension vector. They are originated from the Z3 orbifold projection of the 6-dimensional 
anti-symmetric representation of SU{A) R-symmetry in A/" = 4 SYM: 

r = 7j((^,')^0')7„ for any g G Z3. (22) 

Here 7^ are the regular representation of (7 G Z3 in U{3N) group, the 6-dimensional 
representation of (7 G Z3 in SU (4) R-symmetry group. In particular, 0" corresponds charge- 



2/3 fields in 6 of SU{ 
to a statement in ref. 



) under Rg transformation, while 0" charge- (—2/3) ones. According 
25l |. the 7g-action on correlation functions © is trivial in the planar 
limit. Moreover, gauge invariance of a composite operator ((7j) requires \L — 2A;|/3 to be 
integer, where k is the number of (j)"". Consequently such composite operators have zero 
charge under Rg, and the ADM or spin-chain Hamiltonian ()17|) at the conformal points 
should be the same as the integrable 5*0(6) spin Hamiltonian in = 4 SCYM. 

Explicitly on the conformal fixed line, one has Xa = = A, then the ADM p7j) is 
simplified to 

A ^ 

A ^ 

1=1 
A ^ 

= I^E(2 + ^M+i-2Pm+i)- (23) 
1=1 

Here in the last line, K (a®6) = {a^ -b) Y2i e*®e*. The rank-2 anti-asymmetric representation 
of SU{A) is related to the vector representation of 5*0(6) via an unitary transformation U. 
Notice that {U (g) U)K{U ® U)^ = K, the ADM ^ is indeed the Hamiltonian of the same 
integrable SO (6) spin chain as in ref. jl^ . 
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B. Other J\f = 1 SYM and Wess-Zumino model 

We have shown that, in our present M = 1 model the ADM of (anti-)holomorphic com- 
posite operators corresponds to a Hamiltonian of integrable SU{?)) spin chain. One may 
wonder whether this resuh can be generahzed to other supersymmetric models with global 
SU{M) symmetry, since the integrability condition for a holomorphic spin chain with any 
SU (M) symmetry is quite simple, involving only the identity and permutation operators. 

For a general M = 1 gauge theory with M chiral superfields, our basic observation is 
the following: The contributions of diagrams Fig 2a and 2c are always proportional to the 
identity matrix in SU{M) indices. The quartic scalar vertex in Fig 2b from the gauge 
interaction (the first term of the Lagrangian (j2))) can be effectively treated as an exchange 
of the D-component of the vector superfield, so its contribution is also proportional to the 
identity matrix in SU (M) indices. The non-trivial SU{M) structure in planar one-loop ADM 
can only come from the quartic scalar vertex in Fig. 2b with superpotential interactions 
(the third term in eq. 0). 

Let us first consider a A/" = 1 SU{N) gauge theory, which has M chiral superfields 
transforming as M of a global SU (M) group and as a certain representation TZ under gauge 
group. In order to large limit makes sense for superpotential coupling yet, we restrict 
the representation TZ in which degrees of freedom of chiral superfields are order to A^^ at 
least. In addition to gauge invariance, the superpotential is required to be S'f/(M)-invariant 
and of degree three (for renormalizability) . Therefore the superpotential has to contain an 
S'f/(M)-invariant tensor of rank 3. Consequently the only possibility is M = 3 and the 
superpotential contains the 5'f/(3)-invariant tensor eabc- It is not hard to see that at the the 
planar one-loop level, the ADM of (anti-) holomorphic composite operators contains only the 
identity and permutation operator P and, therefore, can be identified with the Hamiltonian 
of an integrable SU (3) spin chain. 

Now let take the limit in which the 't Hooft coupling tends to zero, then this model 
becomes a Wess-Zumino model with three chiral superfields, each transforming as represen- 
tation TZ of SU{N) group and together as 3 of SU{3) flavor group. With an S'[/(3)-invariant 
superpotential, the ADM of the (anti-) holomorphic composite operators again correspond 
to an integrable SU (3) spin chain, in the same way as the limit we discussed at the end of 
Sec. IV. 
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C. An AT = 2 SYM 



jit complicated. 
The gauge 



If the SYM has an extra U{1) R-symmetry, the story will become a little 
As an example, let us consider the orbifolded M = 2 model proposed in 
group of the model is SU{NY^^ x SU{NY^^ ■ ■ ■ SU{NY^\ The bosonic fields in the vector 
multiplet are denoted as {A^^j, 0/) with I = 1, ■ ■ ■ , K. The matter fields are hypermultiplets 
Bf, (a = 1,2), where Bj belongs to the bi-fundamental representation {N^^~^\ N^^^) of the 
(/ — l)-th and J-th gauge groups, and they form a doublet (labeled by the index a) under 
R-symmetry f/(l)_R x S'f/(2)/j. In order to construct a renormalizable, gauge-invariant and 
S'f/(2)j:j-invariant superpotential, one has to requires K = 2. The superpotential is of the 
form 

W = ^eaMB1B',<^i + B^B',<^2), (24) 

where a, 6 = 1, 2 are SU {2)r indices. 

We can form three types of gauge-invariant, holomorphic composite operators Oi, {i = 
1,2,3) with L sites, which are closed under planar one-loop renormalization: Oi consisting 
of the doublet Bjs only, O2 of (pi only and O3 with mixed Bj and (pj. The closure and 
holomorphy of O2 restricts one-loop planar diagram corrections to O2 to be always diagonal. 
Its anomalous dimension is 

L 



(A„-A,), (25) 



where Xj denotes the 't Hooft coupling of the Jth gauge group. It is interesting to note that 
at conformal points, \j = \h, the operator O2 is protected. 

The Cab in the superpotential forbids the trace operator K to appear in the ADM 
of operators Oi. Moreover, the quartic scalar vertex derived from the superpotential ()24|) 
contains a term (-Bi-B| — BIBI){B2BI — BlBf). The permutation operator P appears, due 
to this term, in the ADM of Oi. 

1 ^ 

= ^ ^(2^^^ - - ^'^^'.^+1)- (26) 
1=1 

Here it should be noticed that (/ > 1) is uniquely determined by Ii due to gauge invariance. 
The ADM in eq. fl26p can again be regarded as the Hamiltonian of an integrable spin chain 
with SU (2) symmetry. 
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The case for is more complicated, since a 0/ field and a neighboring Bj may exchange 
in planar one-loop diagrams. Effectively we can treat an operator as insertions of 0's 
in a Oi operator. Noticing an O3 consists of L Bjs and k 0s is closed with fixed L and k, 
we can write down its ADM explicitly. For example, let us consider a (pi. inserting at i-th 
site of Oi, so that we get operator O3 with L + 1 sites which consists of L Bjs and a 0. Its 
ADM can be written as 

1 A A 

Z=l ! = 1, 

where the permutation operator P' exchange (pi- and neighboring Bj. 

We can see the ADM ()27|1 is similar to the ADM if we suppress the gauge group 
indices. Therefore, we can formally define a Hilbert space n = nr=i ®V with r = © C. 
That is, to put Bj and (pj into a triplet and to consider the set of gauge invariant operators 
{C} = {Ci} © {O2} © {C3}. The ADM for this set of operators is 

1 ^ 

The above ADM can be regarded as the Hamiltonian of an integrable spin chain with SU{3) 
(instead of SU(2)) symmetry. It implies that in the ADM of holomorphic composite scalar 
operators, its symmetry (as spin chain) is enhanced from SU{2) i?-symmetry to SU{3), if 
gauge group indices are suppressed, as allowed by gauge invariance constraint. The difference 
between the Hamiltonian (j^Hj) and ()19|)is only that they have different constant terms. 

Finally, since this M = 2 model is obtained by deforming an Z2 orbifolding of A/" = 4 SYM, 
the ADM for gauge invariant composite scalar operators coincides with that in A/" = 4 SYM 
at conformal points. The same as in the M = 1 case, it corresponds to the Hamiltonian of 
an integrable spin chain with 5*0(6) symmetry. Consequently for the ADM of SYM theories 
we obtain a cascade of integrable structures from orbifolding (or taking quotient of) A/" = 4 
SYM (see Fig. 3). 

D. A conjecture 

In the above examples we have seen that in these orbifolded daughters of TV = 1,2 
SYM, the variation of 't Hooft couplings affects only the total sum of constant terms in 



16 



N= 4 SCYM 




away from conformal line 



N= 1 SYM 



N= 2 SYM 



SU(3) spin chain 



SU(2) spin chain 



FIG. 3: The cascade of integrable structure in SYM. 

the Hamiltonian of the quantum spin chain, and the change by an overall factor in the 
superpotential couplings only leads to an overall factor for the spin chain coupling. Neither of 
them affect the integrability of the spin chain Hamiltonian. So it is the symmetry structure 
of the superpotential that dictates the integrable structure and, moreover, the integrable 
structure should exist for all deformed orbifolded daughters of A/" = 4 SYM with unbroken 
supersymmetries. 

More precisely, we are led to make the following conjecture: In all A/" = 1, 2 orbifolded 
daughters Q, Q] of the AT = 4 SYM, (with a non-abelian global symmetry), with their 't 
Hooft and superpotential couplings deformed away from the conformal points while keeping 
the global symmetry of the superpotential, there is always a nontrivial integrable structure 
in operator mixing at planar one-loop level for (anti-) holomorphic composite operators, at 
least for those consisting of purely scalars without derivatives. Moreover, this integrable 
structure survives in the resulting Wess-Zumino models when all 't Hooft couphngs are sent 
to zero. 



VI. SUMMARY 

We have presented examples in which an integrable structure appears in four dimensional 
A/" = 1, 2 super Yang-Mills theories or Wess-Zumino models, even with conformal symmetry 
broken. In these examples the planar one-loop anomalous dimension matrix for composite 
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operators consisting of (anti-)holomorphic scalars can be written as the Hamiltonian of an 
integrable spin chain with SU{3) (or SU{2)) symmetry. It indicates that the origin of the 
integrable structure in four dimensional Yang-Mills theories is of more profound origin, not 
restricted to superconformal symmetry. Though these spin chains can be formally viewed 
as a subsector of the integrable SO (6) chain found in = 4 SCYmQ], as will be shown in 
appendix B, the operator mixing in our models is completely different, because each scalar 
field factor in the composite operators now carries an extra discrete index A, which is absent 
in A/" = 4 SYM. Moreover, the appearance of an integrable structure in various SYM theories 
enables us to use methods in one-dimensional integrable models to study properties of SYM. 
In appendix A and B below, we will show some examples for how to find one-loop anomalous 
dimensions of composite operators of SYM via solving Bethe ansatz equations. 

Though the study of the ADM of composite operators in A/" = 4 SYM was first moti- 
vated by the BMN limit, which is conjectured to be dual to string theory in the pp-wave 
background, the integrability of the ADM, as a spin chain Hamiltonian, has nothing to do 
with the BMN limit. All our examples confirm this point once more. 

As we have seen both in the A/" = 4 SYM and our models with A/" = 1,2, the length of 
the integrable spin chain can be finite. Moreover we have shown that non-trivial integrable 
structure of these spin chain is actually encoded in the superpotential interactions, rather 
than gauge interactions. In the cases we studied, it is the symmetry of the superpotential 
that determines the symmetry of the integrable spin chain, and it is the superpotential 
coupling appears in the non-trivial part of the integrable Hamiltonian, even though our 
deformation of the superpotential strength has made the theory away from the conformal 
points. 

Recently there are further developments on the A/" = 4 SYM integrable super spin chain: 
A study on the two-loop corrections to the integrable Hamiltonian yielded i ntrig uing evidence 
that the higher order corrections do not break the integrability property [2^. In addition, 
it was shown that there is a relation between the infinite-dimensional non-local symmetry 
of type IIB superstring in AdS^ x and a non-abelian and nonlinear infinite-dimensional 



Yangian algebra for weakly coupled SCYM 



It would be interesting to address similar 



issues in the framework of A/" = 1,2 supersymmetric models, such as those we have presented 
in this paper. Moreover, from the string theory point of view, the following questions arise 
naturally: Is some of the gauge theories studied here dual to gravity or string theory? Could 
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the integrable structure in the supersymmetric gauge theories be originated from or related 
to string theory? 
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APPENDIX A: GENERALIZED HEISENBERG (ANTI-)FERROMAGNET 

We will show that the spin Hamiltonian in Eq. (fT^ is just one of the generalized Heisen- 
berg ferromagnet models discussed in ref. 3l|. For readers' convenience, in this appendix 



we give a brief review on the generalized Heisenberg (anti-)ferromagnet — a quantum sys- 
tem with M components on a one-dimensional chain with nearest-neighbor (short-range) 
interactions. The complete space of state is 

L 

n = (^ni, ni = c''. (ai) 
1=1 

The Lax operator, associated with the n-th site on the chain, L„ a(/i), acts on the tensor 
space Hn ® K with auxiliary space Va = C*^. It is of the form 

LnAl^) = (^il^)In,a + b{fl)Pn,a, (A2) 

where Pn^a is the permutation operator acting on 7i„ ^Va,fi the spectral parameter, and 
a(/i) + = 1, a{fi) = + ie) with e = ±1 corresponding to anti-ferromagnet and 
ferromagnet respectively. 

The transfer matrix TL,a(yu) defined by 

TL,a{fJ') = LlAf^) ■ ■ ■ Ll,a{fl') (A3) 
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is a monodromy around a circle (assuming the periodical boundary conditions with TCl+i = 
Til for the chain). It satisfies the following relations 

where RabifJ') = b{^) + a{fi)Pa^b. Taking the trace on the auxiliary spaces Va and Vb, we get 
the commutative relation 



t(z/)] = 0, t{fi) = tla TlM. (A5) 

This allows us to treat t(/i) as a generating function of commuting conserved quantities: 



M« = z(^-j ln[t(/.)t(0)-i]U=o. (A6) 

By definition, the momentum operator P on a lattice, P = — ilnt(O), is related to the 
discrete shift operator (by one site) 

t{0)ni ®'H2®---®'Hl = 'H2®---®'Hl®'Hi. (A7) 

While the Hamiltonian of the system is 

L 

if = M«=5 5^(PM+i-l). (A8) 

1=1 

Then we can see that for very large L the spin Hamiltonian (jl9p for our M = 1 model is 
nothing but Eq. ()A8|) with e = —1 (for a ferromagnetic systein) 
For Tin = C^, the eigenvalues A(/i) of the operator t(/i) is 



j=i a{fi - /i} i=i ^^^) j=i a{fi -fij')) i=i 



Here the rapidity variables yUj and fi^^^ satisfy the algebraic Bethe Ansatz equations (ABAE): 

n m / (1 



J]a(/if)-/xfc) = n 77771^ [Th 



fc=i k=i (^{fJ'k f^j ) 



n —1^ — = n 'r'~^ V (aio) 
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The eigenvalues of the momentum operator P and Hamiltonian H are, respectively, 



j = l j = l rj 

n n 



Introducing new rapidity variables fiij and fi2,j via 



1 i 

l^? = lf^2,j - le, (A12) 



- 2 

we get the usual expressions for the ABAE (jAlOj) 

_ TT- /ilj - /il,fc - T-T- jUij - /i2,; + 

M /^ij - + 2«£: fj^ yUlj - /i2,/ - 
and for eigenvalues of total momentum and energy, we have 

3=1 ^ ^ 

APPENDIX B: ANOMALOUS DIMENSIONS FROM BETHE ANSATZ 

In this appendix we present several concrete examples to show how to obtain anomalous 
dimensions by solving Bethe ansatz equations. The discussion will be similar to that in refs. 
14 , , but there are new issues to address, associated with the discrete index A arising 
from orbifolding, which is absent in the Af = 4 case. In section 4 we have shown that the 
ADM for holomorphic operators in our A/" = 1 model gives rise to an SU{3) spin chain, which 
can be viewed as a closed sector of the SO (6) spin chain in the Af = 4 SYM. However, as 
we will see below, the operator mixing in our model is completely different, because each 
0-field factor now carries an extra discrete index A that labels the gauge group factors. 
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The SU{3) symmetry of our spin chain is related to Lie algebra A2, which has two simple 
roots, 

'3 1 



«i = U/2'-y|^ «2 = (0,y2), (Bl) 

and the highest weight vectors that generate the fundamental and anti-fundamental repre- 
sentation are, respectively. 

We always take e = — 1 in the rest of this appendix. Then the ABAE ()A13|) follows from 



the general ABAE 



3: 



fiqj - lag ■ w) . _ ^^^^ _ . 5^ 11 11 . _ _ . • 

1. Physics 

Each factor in the composite operator corresponds a site in the spin chain, which is 
occupied by a scalar field, being one component of the SU{'i) triplet 0". For the sake of 
convenience we denote them by (Z, W) in this appendix. 

a ] The ground state 

Because the eigen-energy of the ferromagnet system ()A8|) . namely Fi in eq. (fT^. is 
non-negative, the ground state VL must have zero eigen-energy and total momentum. 
Consequently, must consist of only one component of (Z, y, W). For convenience 
we assume it consists of Z only ^. Since there is no impurities of /ii and /X2, it is not 
hard to see that P^^ = E^^ = 

In our M = 1 SYM, the anomalous dimension of VL is given by the "zero-point energy" 
of the system (Fq in Eq. (HH))): 

ln = Eo = -^(^1 + A2 + A3 - 3A,). (B4) 

Along the conformal line, = Xa, one has 7^^ = 0, so that the dimension of Q is 
indeed protected by super conformal invariance. 



[2] In Af — 4: SYM this state correpsonds to the ground state of the BMN operator. 
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b ] States with impurities 

In the present case there are two types of impurities in the spin chain, labeled by 
two rapidities: /iij and fi2,j, which are associated with the two simple roots di and c52 
respectively. The states with impurities correspond to excitations of W and Y (i.e. the 
replacements of some Z hj W and/or Y) in the ground state Q. Since now both Wi — di 
and Wi — d2 are also weights, in addition to purely fii or /i2 impurities, /xi— /i2 bounded 
impurities are also allowed. However, the weight wi — 02 is not equivalent to wi, so it 
does not lie in the fundamental representation. If we restrict ourselves to holomorphic 
operators, a single /i2-impurity without being bounded to a /xi-impurity on the same 
site is not allowed. The physical interpretation is the following: A single /xi-impurity 
(wi — di) creates a W replacement in the state Q, while a yUi — /i2 bounded impurity 
(wi — di — d2) creates a Y replacement. But an individual /i2-inipurity (wi — a2) 
would create a Y replacement in Q, breaking the holomorphic nature of the composite 
operator. 

c ] The trace condition 

An important observation is that the composite operators that we are considering are 
gauge invariant after taking the trace in the product gauge group U{N) x U{N) x U (N), 
while the scalar fields belong to bi-fundamental representations of two gauge groups, 
instead of the adjoint of one group. So gauge invariance and holomorphy of the 
composite operators requires that we are dealing with a chain with length L = 3k 
with k integer, and 



In other words, after shifting the chain by one site three times in the same direction, 
we should obtain the same composite operator. Therefore we have the cubic trace 
condition: 



It is easy to verify that this trace condition is consistent with the ABAE (IA13|) . It 
follows that the total momentum can be only 2n7r/3 with n integer. 



(B5) 




(B6) 
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We note that the condition ()B5|) or ()B6|) corresponds to, but is very different from, 
the trace condition of the 5*0(6) spin chain in A/" = 4 SYM Q|. In the latter case, 
one has power 1 instead of power 3 in Eq. ()B5|) and Eq. ()B6|) . The reason for power 
3 in our model is obviously related to the orbifolding by Z^, which leads to an extra 
three-valued index A for the holomorphic scalars. Below we will see that it is the 
appearance of this index, though the spin Hamiltonian is "blind" to it, that makes the 
operator mixing very different from that in A/" = 4 SYM. 



2. One impurity 

The first non-trivial and interesting case is a single /ii-impurity in Q. This is the simplest 
example which shows how and why the operator mixing in our model is very different from 
that in the Af = 4 SYM, though the spin Hamiltonian may be viewed as a closed subsector 
in the SO (6) spin chain in the latter. In A/" = 4 SYM, a single impurity with non-zero 
momentum is not allowed, because of the trace condition. As shown below, however, the 
present case does allow a single /ii-impurity. 

The trace condition and the ABAE now reduce to 

V/^i - V 

Here we have used the fact L = 3k with k integer. The above equation yields 

2mT 



E = e(/ii) =4sin2^. (B8) 



Then anomalous dimension is 



In the language of A/" = 1 SYM, we are now considering excitations with one W replace- 
ment in the ground state Q. There are three distinct possibilities: 

O, = ii{WiZ2Z^{ZiZ2Z;f-^}, 

02 = tT{Z^W2Z^{ZiZ2Z;f-^}, (BIO) 

03 = ti{Z^Z2W^{Z^Z2Z^f-'}, 
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where the subscripts 1, 2, 3 are values of the index A, labehng three types of bi- 
fundamentals. In the above operator basis, the ADM ()19|) reads 

/ 



87r2 



M + To, 



M 



\ 



2 -1 -A 
-1 2 -1 
-1 -1 2 



(Bll) 



The matrix M has eigenvalues {3, 3, 0}. Therefore we obtain the eigenvalues of F as follows: 



7i = 72 



73 = 



(B12) 



The results are precisely the same as one obtains from the ABAE. The corresponding eigen- 
vectors are {aOi,b02, —{a + 1)0^} for 71 = 72 and c{0i,02,0^} for 73 with arbitrary 



constants a, 6, c. 



3. Two impurities 

Two impurities can be either two /ii or one /ii and one ^2- Let us consider two yUi- 
impurities first, i.e. excitations with two W replacements in the ground state. The trace 
condition and the ABAE now reduce to 

/^l.l + ^ . /^1.2 + i ^ ^2inn/3 

.Ati,i-V - Ati,2 - 

For n = 0, we must impose /ii^i = — 2, then only real solutions are allowed. We get the 
momenta = 2mTX / (L — 1) and the anomalous dimensions (with m an integer): 

7r°) = ^sin^f^ + 7,. (B14) 
For ri = 1, 2, however, both of real and complex solutions are allowed. If yUi 1 is real, from 



Eq. (jB13j) we obtain 

(L - 1)^(^1,0 = mTT + T /ii,2 = (B15) 

v3 V3/Ui,i=f1 

where = arctan(x), the integer m parameterizes different branches of the logarithmic 
function. Eq. ()B15|) in general can not be solved analytically. For very large L, however, 
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we have approximately 9 = -^{fii^i) ~ mn/L. Consequently the anomalous dimensions are 
given by 



7^^'^) = ^(l + ^cos2^T^sin2^)+7, 



\h , 1 2m7r 2rmx ^ 1 s n 

= ^(1 + i ^ T ^ ^) + 7. + O(^). (B16) 

Now let us consider complex solutions. Notice that for L ^ oo, the LHS of the second 
equation in Eq. ()B13|) grows (or decreases) exponentially if Im/ii i 7^ 0. Hence the RHS of 
this equation together with the first equation in Eq. ()B13|1 lead to the solutions 

= + or = (B17) 

and the anomalous dimension 

0^^i^ + 0„. (B18) 

Notice that 7 < min(7m~^'^''). It indicates that complex solutions correspond to bound 
states. 

Next we consider the bounded impurity of one /xi and one /i2- The ABAE together with 
trace condition now reduce to Eq. ()B7|1 plus 

i^l^J^ = 1. (B19) 

The solution from those equations are /ii = cot ^ and /i2 = 00. It yields the following 
anomalous dimension 



= FT^(^i) + 7n = TTT ^ + 7o- (B20) 



The result is the same as in Eq. ()B9jl . It just reflects the fact that replacing a Z in the 
ground state either by W or hj Y leads to the same anomalous dimension. 

4. The highest excited state 

For a finite chain the excited state with the highest energy contains as many as possible 
impurities. Taking the logarithm of the ABAE ()A13|) we have 

^711 n,2 
k=/=j k=l 
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feT^j fc = l 

where we have used that fact that the discreteness of the Bethe roots requires them to be 
pushed to different branches of the logarithm function. In general the ABAE ()B21|) can 
not be solved analytically with more than two impurities. In the thermodynamical limit 
L — > oo, however, the Bethe ansatz equations are simplified significantly j2J|. In this limit 
j/L is replaced by a continuous variable x, and the ABAE ()B2H) is replaced by a set of 
integral equations: 

= nx + jdy .Jdy i}{^,{x) - /.i(y)). (B22) 

Taking derivatives with respect to fii and fi2, we have 

2 , -I = vrpi(/ii)+ / d/i — , ^ - / rf/i 



+ l """" i-oo (pi-p')' + 4 y_oo + l 



= 7rp2 p2 + / rf/i 7 -K^mi ~ / 7 /N2 I 1 ' B23 



where the densities pi and p2 are defined by 



Pi(/^i) = T^TT' ^2(^2) = -r~T^- 



The equations ()B23|) can be solved by means of Fourier transformation. The results are as 
follows: 

dk jj^,^. 2 cosh k 



This yields 



271 4coshH-l' 



2 [°° 1 

dx pi(x) = -, / rfx p2(x) = -. (B26) 



Consequently we have 2L/3 pi-impurities and L/3 p2-impurities for the highest excited 
state. They fill all sites on the ferromagnet chain. In our M = 1 SYM, it implies that there 
are equal number of W , Y and Z scalar fields in the composite operator. Recalling L = 3k 
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with k integer and denoting the SU{3) triplet as (p"", (a = 1, 2, 3) again, the operator has 
the following form: 

O ~ tr(ea6,0>V')'- (B27) 

It is an SU (3) singlet with zero total momentum. 
The anomalous dimension of this operator is 

It grows linearly with L, the length of the chain. 



Similar to the procedure in 
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we can also calculate the anomalous dimensions of 



operators with a few W and Y replacements. But let us stop here. 
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